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Abstract. We present the detection of molecular oxygen with Odin toward the dense molecular 
core pOphA, which is part of a region of active star formation. The observed spectral line is 
the (Nj = li — lrj) ground state transition of O2 at 119 GHz (2.5 mm wavelength). The center 
of the line is at the LSR velocity of a number of optically thin lines from other species in the 
region and the O2 line also has a very similar, narrow, line width. Within the 10' beam, the 
line intensity is J TAdv = 28 mK kms , which corresponds to 5<r of the rms noise. A standard 
LTE analysis results in an O2 abundance of 5 x 10~ 8 , with an uncertainty of at least a factor of 
two. We show that standard methods, however, do not apply in this the coupling of the 

Odin beam to the source structure needs to be accounted for. Preliminary model results indicate 
O2 abundances to be higher by one order of magnitude than suggested by the standard case. 
This model predicts the 487 GHz line of O2 to be easily detectable by the future Herschel-HIFI 
facility, but to be out of reach for observations on a shorter time scale with the Odin space 
observatory. 
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1. Introduction 

The energy balance of the interstellar medium (ISM) relies on a complex interaction 
of various heating and cooling processes, which of course also affect the chemistry being 
activated in the medium. This in turn regulates the amount of the major cooling species 
in the gas. Heating procsses include mechanical and radiativ energy input, whereas, 
in dense clouds, the cooling is dominated by molecules and dust particles, resulting in 
low equilibrium temperatures on short time scales. Whereas the coolant CO is readily 
observable from the ground, water and oxygen molecules need to be observed from space. 

Odin is thus a mission in orbit around the Earth dedicated for the observation of H2O 
and O2 and the facility is shared by aeronomers and astronomers on an equal time basis. 
The millimeter/submillimeter (mm/submm) observatory Odin is described by Frisk et 
al. (2003) and Olberg et al. (2003), whereas Odin-astronomy is reviewed by Hjalmarsson 
et al. (2003). These articles are all collected in the special Odin edition of Astronomy & 
Astrophysics, which of course also contains numerous papers on initial results obtained 
with Odin. Energy level and excitation diagrams for the principal observable molecules 
can be found in the review by Liseau (2001). 

Prior to the launch of Odin in February 2001, it was already clear from the observations 

f Odin is a Swedish-led satellite project funded jointly by the Swedish National Space Board 
(SNSB), the Canadian Space Agency (CSA), the National Technology Agency of Finland (Tekes) 
and Centre National d'Etude Spatiale (CNES). 
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by the Submillimctcr Wave Astronomy Satellite (SWAS) that H2O was significantly less 
abundant in the ISM than was previously thought. Even more intriguing, perhaps, was 
the apparently total absence of O2 (Goldsmith et al. 2000). However, a couple of years 
later, the tentative detection of the 487 GHz line (see Fig.^) toward the dense molecular 
core pOphA was announced by Goldsmith et al. (2002). 

The SWAS result was challenged by Pagani et al. (2003), who were unable to confirm 
the presence of oxygen in p Oph A based on observations with Odin in the ground state 
transition of O2 at 119 GHz ( Fig.^J. Taking into account the fact that this involved 
a different line at considerably lower angular resolution toward a somewhat different 
position, rendered the status of the SWAS result inconclusive. 

Besides the negative result for the dense cloud pOphA, Pagani et al. (2003) reported 
also on unsuccessful observations of nearly a dozen galactic sources in the O2 (Nj = 
li — lo) line at 119 GHz. Furthermore, Wilson et al. (2205) recently published an upper 
limit in this line toward an extragalactic object, viz. the Small Magellanic Cloud (SMC). 
Contrasting with these negative results, we report here the detection with Odin, at the 
5fi level, of the 119 GHz line of O2 toward the dense cloud pOphA. 

2. Odin O2 Observations 

The Odin observations of p Oph A were performed during 2002 and 2003 and the total 
on-source integration time was 78 hours. At 119 GHz, the half-power-beam-width of the 
1.1m telescope is 10'. As back-end, a digital autocorrelator of 100 MHz bandwidth and 
125 kHz resolution was used. In Dopier- velocity units, this corresponds to 250kms _1 
and 0.3kms _1 , respectively. More exhaustive details regarding these observations and 
the subsequent complex data reductions are provided by Larsson et al. (2005). 

In Fig.0 showing the final stage of this data reduction, the Odin 119 GHz spectrum is 
presented. At full spectral resolution, the O2 line occupies five autocorrelator channels, 
which in the figure have been binned into one channel. This is equivalent to the integrated 
line intensity, §T^dv z = (28 ± 5.5)mKkms _1 . 

Other Gauss-fit parameters are: peak Ta — (17.5 ± 3.5) mK, line center velocity 
v lsr = (3-0 ± 1.0) kms -1 and line width A^fwhm = (1-53 ± 0.03) kms -1 . The ve- 
locity and width of the O2 line are in agreement with that of the PDR recombination 
line emission observed by Pankonin & Walmsley (1978) with a 7'.8 beam, i.e. comparable 
to the spatial resolution of Odin. As also indicated by the C 18 observations of Frerking 
et al. (1989) with a 2'3 beam, the line parameters center velocity and width of optically 
thin species appears constant on various angular scales. The values of both parameters 
are significantly different for the 487 GHz line of Goldsmith et al. (2002). 

The Odin observations at full spectral resolution, and the comparison with other melec- 
ular line data, will be presented by Larsson et al. (2005). 

3. O2 Abundance 

Of general interest is the relative concentration of various forms of oxygen in the inter- 
stellar medium and in particular that of the molecular constituents. Below, we shall ex- 
amine the derivation of the molecular abundance in some detail, as widely used standard 
methods will encounter difficulties in this particular case, leading to erroneous results. 

3.1. Column density of oxygen 

The data for radiative (A u i) and collisional (q u i) de-excitation have been provided by 
Marechal et al. (1997) and by Bergman (1995 and private communication). The collisional 
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Figure 1. Rotational energy level diagram of O2 for upper energies below 30 K. Transitions 
which are accessible to various missions discussed in the text are also indicated (adopted from 
Liseau 2001). 



excitation Nj = 1q to N = J, such as (lo — li), is forbidden and we assume that the 
collisional population of the li level occurs mainly via 32 and hence, the critical density 
for 2 (Nj = li - 1 ), viz. 



n-crit = An, 10/932, 11 (?k) 



(3.1) 
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Figure 2. The 5 channels of the O2 feature have been binned into 1 channel, corresponding 
to the integrated line intensity jT\dv at the 5 a level. All other features in the entire useful 
spectral range over 180 km s -1 are below the 3 a leve, shown as the dashed horizontal line. 

becomes > 10 3 cm" 3 but < 10 4 cm~ 3 for T k < 100 K and is likely lower than densities 
in the cloud. Local Thermodynamic Equilibrium (LTE) provides thus a reasonable as- 
sumption for the population of the lower levels. The weak O2 (Nj = li — lo) line is most 
likely optically thin and can be expected to exhibit a nearly Gaussian shape, reflecting 
small scale stochastic motions in the cloud. 

The observed integrated intensity of the optically thin 119 GHz line, formed in LTE, 
is then given by 
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where /b is the beam filling, 7y m b is the main beam efficiency and Tr is the radiation 
temperature of the source. In LTE, the only dependence is on the temperature, which for 
our present purposes is well represented by a power law, i.e. for the 119 GHz transition, 
the function $(7k) is given by 
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in units of Kcm 3 s _1 . In the above equations, the temperature of the back-ground 
radiation field is Tb g = 2.734 K, the statistical weight of the upper level is gu = 2J+ 1 = 
3, the temperature of the transition is Tu t io — 5.70 K, the quasi-Planck function is 
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Figure 3. Continuum observations of the p Oph cloud at 200 pm with the 2' beam of the bal- 
loon-borne submm-telescope PRONAOS (Ristorcelli et al. 2005). pOphA is the hot spot above 
the center of the figure. Intensities, I v , axe given in MJysr -1 (see the scale bar to the right). In 
addition to the 200 /im-map shown here, full 40' x 40' maps have also been obtained at 260, 360 
and 580 pm. 

Ji/{T) — Tn j io/[exp(Tn i io) — 1] and the partition function is approximated by Q(Tk) ~ 
kT k /hB, where the rotational constant is B = 43100.460 MHz. 

Expressing the line intensity in Kkms -1 , the column density of O2 is then simply 
found from 

N{0 2 ) 10 15 T k - 67 (cm" 2 ) (3.4) 

On the scale of the Odin beam, temperatures have variously been estimated from gas 
tracers (e.g., CO; Loren et al. 1983) and multi-wavelength measurements of the dust 
emission (e.g., Ristorcelli et al. 2005), yielding an average of about 30 K and, hence, 
N(C>2) ^ 1 x 10 15 cm -2 . As can be seen from Eq. 3.4, the column density of O2 is not 
very sensitive to the temperature and, e.g., an uncertainty of ±10 K would amount to a 
change in A^(0 2 ) by no more than 20%. The critical parameter for the determination of 
X(0 2 ) obviously is jV(H 2 ). 

3.2. Column density of hydrogen 

To find the oxygen abundance, X(Oa) = N(02)/N(H.2), we need to specify the column 
density of hydrogen, H 2 (see: Eq. 3.2). Based on the C 18 maps of Tachihara et al. 
(2000), Pagani et al. (2003) estimated for Odin O2 observations with a 10' beam the 
H2-column to be A^Ii^) = 4 x 10 22 cm~ 2 . This is higher than values recently obtained 
from [CI] observations by Kulesa et al. (2005) on a much smaller, i.e. 3'5, scale, viz. 
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Figure 4. Theoretical (LVG) O2 spectrum between 100 ^tm and 1cm for pOphA, based on the 
Odin observations in the O2 (Nj = li — lo) line of p Oph A. The uniform source is assumed to fill 
the 10' beam. The bar shows the frequency coverage of Herschel-HIFI, admitting the 487 GHz 
(33 — I2) line at its low-frequency end (see also Fig.0, whereas the stronger 425 GHz (32 — I2) 
line falls unfortunately outside the HIFI band. 



(2 — 3) x 10 22 cm -2 and those derived from submm dust continumm observations by 
Ristorcelli et al. (2005), viz. 1 x 10 22 cm' 2 for the 10' Odin beam. The earlier map of the 
LTE column density of C 18 by Wilking & Lada (1983) would indicate a one order of 
magnitude higher N(H2). This is also the value obtained, on smaller angular scales, by 
Kamegai et al. (2003) from a large [CI] map of the p Oph cloud. 

Various assumptions come into these estimates, such as the abundance of C ls O or the 
opacity of the dust grains, leading to a considerale spread in 7V(H2). For the moment, 
we shall adopt a weighted average of 2 x 10 22 cm -2 , a value which is associated with an 
uncertainty of at least a factor of two. 
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3.3. 02 Abundance: extended source 

Applying the estimates of temperature and H2-column density to the O2 observations 
with Odin yield an oxygen abundance X(02) ~ 5 x 10~ 8 , if we assume that both rj m b 
and /b are close to unity. The corresponding rotational O2 spectrum between 100 firm 
and 1 cm is displayed in Fig.|3| From that we infer that Odin is not capable to detect 
from pOphA the other, in principle accessible, O2 line, i.e. the (33 — I2) line at 487 GHz 
(2f3 beam). It is also clear that the detection of the 487 GHz line with Herschel-HIFIt 
(45" beam) would require the O2 source to be significantly more compact than 10'. 

3.4. O2 Abundance: compact source 

The value of X{0 2 ) derived above is based on an analytical model of the line excitation, 
where one has made the implicit assumption that the emitting region is an isothermal and 
homogeneous plane-parallel slab of gas. These idealized conditions would have to apply 
on the angular scale of 10', corresponding to a physical length scale of half a parsec at the 
distance of the p Oph A cloud. This picture is inconsistent, however, with observations on 
smaller scales of the gas and the dust (see, e.g., Ashby et al. 2000 and Johnstone et al. 
2000, respectively), which show pOphA to have a centrally condensed core structure. In 
other words, strong gradients in temperature, density and velocity field are known to be 
present inside the Odin O2 beam. In addition, UV-radiation illuminates and heats the 
core preferentially only from one side (Liseau et al. 1999). 

Modeling the complex structure of this core is in progress. Preliminary results suggest 
that the line formation region of O2 119 GHz is relatively compact, with a size of 22". 
Consequently, the implied abundance is relatively high, A(02) = 6 x 10~ 7 , and the 
observed weakness of the line is due to heavy beam dilution (roughly 0f ource /^ cam = 
10~ 3 ). Predictions of this model for the 02 487GHz line are for Odin (2'3 beam) a 
peak Ta=22uiK and for Herschel-HIFI (45" beam) 0.2 K. Whereas the observation of 
the former would imply prohibitive integration times, the latter would be a matter of 
minutes. We look impatiently forward to the launch of Herschel in the 2007/08 time 
frame. 



4. Possible implications 

The weak detection of O2 by Odin will likely not add to the present discussion of 
the energy balance of the ISM, as the role of a coolant of importance had already been 
altered before. For instance, the total cooling in the O2 lines of the spectrum of Fig.U 
amounts to a mere 7 x 10~ 5 Lq. 

It is conceivable, though, that this detection eventually will lead to an increased un- 
derstanding of the O2 chemistry of the dense and cold ISM and that chemistry models 
may benefit from the Odin achievement. 

5. Conclusions 

In summary, we briefly conclude the following: 

i. Very deep integrations (3.5 mK rms) with Odin in the 1119 GHz ground state line 
of O2 toward the dense cloud core pOph A resulted in a 5a detection of the line. 

ii. Standard analysis of an optically thin line formed in LTE leads to the beam-averaged 
column density of molecular oxygen for pOphA, viz. A(02) = 1 x 10 15 cm -2 . 

f http: //www. sron.nl/divisions/lea/hifi/ 
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iii. Estimates of the oxygen abundance need to rely on assumptions regarding the 
source size and structure with respect to the 10' beam of Odin at 119 GHz. This results 
in X(0 2 ) = 5 x 10~ 8 for a homogeneous extended source and X(0 2 ) = 6x 10~ 7 for a 
more realistic source structure, respectively. 

iv. This issue will most likely become resolved by observations with the future Herschel- 
HIFI, for which we predict an O2 intensity of a few tenths of a Kelvin at 487 GHz. 
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